Atom interferometric detection of the pairing order parameter in a Fermi gas 
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We propose two interferometric schemes to experimentally detect the onset of pair condensation 
in a two spin-component Fermi gas. Two atomic wave-packets are coherently extracted from the gas 
at different positions and are mixed by a matter-wave beam splitter: we show that the spatial long 
range order of the atomic pairs in the gas then reflects in the atom counting statistics in the output 
channels of the beam splitter. Alternatively, the same long range order is also shown to create a 
matter-wave grating in the overlapping region of the two extracted wave-packets, grating that can 
be revealed by a light scattering experiment. 
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The experimental possibility of controlling at will the 
scattering length a between two spin components of 
fcrmionic atoms via a Feshbach resonance has opened 
the way to a comprehensive study of the pairing transi- 
tion in a degenerate Fermi gas [By Q The weakly inter- 
acting limits are well understood theoretically: the phase 
transition is the Bose-Einstein condensation (BEC) of di- 
atomic molecules (a = + ) or the BCS transition due to 
pairing in momentum space (a = 0~). But one can now 
investigate experimentally the theoretically challenging 
crossover region, including the unitary limit \a\ = oc 

While the standard techniques used for atomic BECs 
have allowed to detect and characterize a molecular 
BEC , a debate is still in progress about experimental 
signatures of pair condensation for a negative scatter- 
ing length. Several proposals have been put forward ||; 
none of them was proved to demonstrate the existence 
of long range order in the pairing parameter. First ex- 
perimental evidences of a condensation of fermionic pairs 
in the crossover regime have been recently presented , 
based on a fast ramping of the magnetic field to con- 
vert pairs on the a < side into bound molecules on the 
a > side, and on the observation of the Bose condensed 
fraction of the resulting gas of dimers. This method is 
expected to work only when the fcrmionic pairs are small 
enough, that is in the vicinity of the Feshbach resonance, 
fejr|a.| > 1 where Uf is the Fermi momentum. 

In this paper, we propose a more direct and general 
way of proving the condensation of pairs, by a measure- 
ment of the pairing order parameter, which is not re- 
stricted to the small pair regime fej?|a| > 1. This pro- 
posal is the fermionic analog of the atom interferometric 
measurement of the first order coherence function G^ 1 ' 
of a Bose gas Q. More subtle schemes than the obser- 
vation of the mean atomic density have however to be 
introduced as there is no long range first order coherence 
for fermions. Their experimental implementation would 
constitute a remarkable transposition of quantum optics 
techniques to a fermionic matter field. 

In the current theories of the superfluid state in 



fermionic systems 0, 0, Hj] , the onset of pair condensa- 
tion is defined by a non-zero long-distance limit x ab = 
\xa — xb| — > +oo of the pair coherence function 

G« r (x A ,x B ) = (^j(x A )#{(x A )* i (x B )^ T (x B )), 

(1) 

this function then factorizing in the product of the order 
parameter in and the complex conjugate of the order 
parameter in x^. In the following, we shall propose two 
distinct methods to measure G^ ir (xa , xg ) , both relying 
on the coherent extraction of two atomic wave-packets in 
tca.b and their subsequent beating. The first method is 
based on a two-atom interferometric technique inspired 
by two-photon techniques |9|: it relies on atom counting 
in the two output channels of a matter-wave beam split- 
ter. The second method is based on the coherence prop- 
erties of light elastically scattered off the matter-wave 
interference pattern of the two overlapping wave-packets. 

Consider a gas of spin-1/2 fermionic atoms at thermal 
equilibrium in a trap. At the time t = 0, the trap poten- 
tial is suddenly switched off and the atom-atom interac- 
tions brought to a negligible strength, so that the subse- 
quent propagation is the one of a free atomic field. At 
the same time, a suitable short pulse of spin-independent 
optical potential is applied (Fig^l to the atoms situated 
in regions of size £ u around the points x^ and x^ so to 
impart them a momentum kick of respectively ko ± ki by 
means of Bragg processes and to produce wave packets 
which are a coherent copy of the field in the trap, but for 
a shift in momentum space, ko is taken orthogonal to 
X/i -xb, while ki is parallel to it. The magnitude Kk\ of 
the counter-propagating momentum kicks is taken larger 
than the momentum width Ap of the gas, which is on 
the order of the Fermi wavevector kp in the resonance 
region (|a| = +00) and in the weakly interacting BCS 
regime (a < 0), or on the order of H/a in the case of a 
molecular condensate (a > 0). The size of the extraction 
region t u is taken much smaller than the distance xab 
between the extraction points. This latter is taken as 
macroscopic, that is much larger than any other length 
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scale of the problem, e.g. the Fermi distance 1/kp and 
the Cooper-pair size £bcs- 

In Heisenberg picture the field operator at the end of 
the optical pulse can be related to the initial one by ^3 : 



momentum ko ± ki will be again spatially separated and 
centered at X± = X + fo(ko ± ki)(i2 —ti)/m. The total 
field can be written as ^ a (x, t 2 ) = (x, t 2 )+$~ (x, t 2 ) + 
$ g (x, t 2 ) where the contribution of each packet is 



* CT (x, At) = it(x - Xj4 ) e^ ko+kl )-( x - XA )* CT (x)+ 

+ u(x - x B ) e^ko-kO-Cx-xi,) ^ (x) + £bg( x) (2) 

The atoms which are left in their original momentum 
state as well as the ones having received a different mo- 
mentum kick during the extraction process are included 
in the background field ^ a s : as they spatially separate 
during the evolution, they will be omitted in the discus- 
sion [13j • We shall assume for simplicity that the extrac- 



tion function is a Gaussian, u(£) 
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of 



size l u \ its peak amplitude uo is of modulus less than one. 

This out-coupling scheme produces two atomic wave- 
packets traveling with momentum ko ± ki starting from 
respectively ~x.a,b- At a time t% = wixab /2ft|kx|, the two 
wave-packets superimpose around X. As mentioned in 
the introduction, the mean density profile in the overlap 
region does not show fringes so that more elaborate ma- 
nipulations have to be performed onto the atoms in order 
to measure the pair coherence function C?p^ ir (x^, x^). 

Atom- number correlations: At t — ti, the two 
overlapping wave-packets of momentum respectively ko± 
ki can be coherently mixed by a spin-insensitive 50-50 
matter-wave beam splitter, with reflection and transmis- 
sion amplitudes of momentum-independent phase differ- 
ence (j>. Such a beam splitter may be realizedyij by ap- 
plying a pulse of sinusoidal optical potential U(x,t) = 
4hfl(t) sin 2 (ki • x + 0/2) Q. At a time t 2 after the 




FIG. 1: First proposed set-up: atoms are extracted by a Bragg 
process from the gas at points xa,b, using pairs of laser beams 
of wavevectors — ko, ki and — ko, — ki respectively; at their 
overlap position X, the two atomic wave-packets are coher- 
ently mixed by a laser standing wave acting as a 50-50 beam- 
splitter with adjustable phase shift (j>\ the number of atoms 
in each wave-packet is measured at the final positions X±. 



splitting procedure, the two emerging wave-packets of 



**(€ + X±,t a ) = 
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*a(x B ,A + • (3) 



1Z(£, t) is the free-particle propagator and 9 is an irrel- 
evant propagation phase which depends on the details of 
the beam splitting procedure. The unitarity of 1Z ensures 
that the results to come do not depend on t 2 . 

The operator N^ 1 giving the number of atoms with 
spin a in the wave-packet ± is obtained by integration 
of ^(x, t 2 )^> a (y:, t 2 ) over the spatial extension of the 
packet ± at time t 2 . The operator giving the atom 
number difference between the two wave-packets is then 
D a = A+ — N~. Its expectation value (D a ) involves the 
first-order coherence function (^(x^ + £) $ ff (x^ + £)) 
of the initially trapped atoms, and therefore vanishes for 
a macroscopic distance xab 3> h/Ap: we shall now take 
= (N~) = N a . Information on the pair coherence 
function G^ h . is obtained from the correlation between 
the two spin components: 
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(4) 



From an experimental measurement of the dependence 
of Cjj, it is therefore possible to determine whether the 
system has long range order or not. 

An explicit calculation of Cfj as a function of the en- 
ergy gap A can be performed by using the zero temper- 
ature BCS theory 0, , with predictions that are accu- 
rate in the weakly interacting limit only. In the large 
xab limit, only the cj> dependent part of Eq.Q) has a 
non-zero value: C-j-j, = cj < |' ) cos(2</>). In the local den- 
sity approximation, and assuming for simplicity that the 
mean densities are the same in the two extraction points 
and in the two spin states, we find an analytical expres- 
sion for a wide extraction region £ u 3> ^bcs ; both in the 
weakly interacting BCS regime 
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and in the regime of a molecular condensate 
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Cj j is obtained by averaging over many realizations of the 
whole experimental procedure starting from a trapped 
gas in the same initial conditions, so that a knowledge 
of the signal-to- noise ratio is relevant. We estimate the 
noise by the standard deviation of D a : from Wick's the- 
orem and to leading order in uo, (-D 2 -) — 2N„ which 
shows that the shot-noise ^l) in the initial extraction 
process is the dominant source of noise. The number 
of realizations over which to average therefore scales as 
{N a /C { ^) 2 , which is on the order of 1 in the BEC limit 
and on the order of (Ep/A) 2 in the BCS limit. 

Light scattering off the matter-wave grating: 
Information on the pairing coherence function G^ ir of 
the trapped gas can also be obtained by means of 
light-scattering off the matter-wave interference pattern 
formed by the overlapping wave-packets at t = t\ around 
point X, which is taken in what follows as the origin of 
the coordinates. As already mentioned, the mean density 
does not show fringes. On the other hand, fringes appear 
in the opposite spin density-density correlation function 

^(X.X') - (*}(x,t 1 )*{(x',t 1 )*;(x',ti)* T (x,t 1 )). 

As a guideline, one performs an explicit calculation for 
BCS theory: one finds that, as soon as a condensate of 
pairs is present, fringes appear as a function of the center 
of mass coordinates (x+x')/2 of a pair, with a sinusoidal 
oscillation of wavevector 4ki . Their amplitude is propor- 
tional to the product of the in-trap anomalous averages 
in x A and xb and extends up to relative distances |x — x'| 
of the order of the Cooper-pair size £bcs ■ This matter- 
wave grating is not easily detected in position space since 
its spatial period is smaller than the mean interatomic 
distance. We therefore switch to Fourier space: 

^f(q,q') = J dxdx'e-^+^'^jf^x'). (7) 

Since fringes show up on the center of mass coordinate 
with a wavevector 4ki we limit ourselves to the region 
q' = q ~ — 2ki [lif. Taking into account the free expan- 
sion during t\, one then obtains: 

G\f (q, q) = e ih A ^' m J d£ d£ e^ta+^O ■(€+€') 
x «*(£)«*(€')«(€ -bM*'-b)x 
( *}(x A +€) *I(x A +0 *i(x B +€'-b) * T (x B +£-b)) 

(8) 

where b = ft(q+2ki)fi/m and Aq — |q+2ki|. Remark- 
ably, for q = — 2ki one recovers the factor in front of e 2 ^ 
in Eq.Q. 

This Fourier component of is detectable in the 
angular patterns of the elastic light scattering from the 
atomic cloud. The incoming laser intensity has to be 
weak enough to avoid saturation of the atomic transition. 
Optical pumping processes have to be negligible during 



the whole measurement time: the mean number of scat- 
tered photons per atom has to be much less than one 
not to wash out the information on the internal atomic 
state. The imaging sequence is assumed to take place 
in a short time so that the positions of the atoms can 
be safely considered as fixed. For each realization of the 
whole experiment, a different distribution of the atomic 
positions is obtained, and consequently a different angu- 
lar pattern for the elastic scattering. Information on the 
density-density correlation function will be obtained by 
taking the average over many different realizations. 




FIG. 2: Second proposed set-up: atoms are extracted from 
the cloud at xa,b and create a matter-wave grating when the 
wave-packets overlap at X; one detects this grating by shin- 
ing a pair of counter-propagating, a± polarized laser beams on 
the overlap region and by beating the two resulting backscat- 
tered light beams on a beam splitter: as function of the mixing 
phase m ix, the beating intensity 7(<^ m ix) averaged over many 
realizations presents fringes revealing the pair long range or- 
der. Dashed (dot dashed) lines: a+ (<r_) polarization. 

We consider here the simple case when the laser field 
frequency is close to resonance with the transition from 
the F g = 1/2 ground state to a F e = 1/2 hyperfine com- 
ponent of the excited state. In this case, a± polarized 
light interacts only with atoms respectively in the |,| 
spin state. The geometry adapted to get information on 
the condensation of pairs is shown in Fig|2 a pair of 
mutually coherent laser beams with a common intensity 
i; nc is sent on the atomic cloud with opposite circular po- 
larizations a± and opposite wavevectors ±ki. We look 
at the mutual coherence of two back-scattered beams in 
opposite directions ±k sc [l6|. with opposite circular po- 
larizations. Within the Born approximation (valid if the 
cloud is optically dilute and optically thin), the ampli- 
tudes in Fourier space of the scattered light on the cir- 
cular polarizations cr± are related to the ones of the in- 
coming field by 

£5(±k sc ) = Ap u (±q,t 1 )E™(±k 1 ), (9) 
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where A is a factor depending on the dipole moment 
of the transition and on the atom-laser detuning, and 
p CT (q, t\) is the Fourier component at the transferred 
wavevector q = k sc — hi of the density operator 
>frj-(x, ti)3 r tr (x, ti) at time t\. This mutual coherence is 
quantified by the correlation function 

I-+/hno = ( [E^(-k sc )Y E s ^(k sc ))/L mc 

= ^([^(-q^i)] 1 ^(q,ti)> = |A| 2 0W(q,q). (10) 

By using (jHJ, one indeed sees that the correlation func- 
tion / |_ can reveal the pair long range order. In the 

weakly interacting BCS regime as well as in the one of 
a molecular condensate, it has a simple expression for 
q~-2k i: 

As a function of Aq, it has a narrow peak with a height 
proportional to the correlation function C-\i of the first 
proposed set-up and with a width l/£] such that l\ — 
1"^ + (hti/m£ M ) 2 . Experimentally, this can be determined 
by beating the two scattered beams: as a function of the 
mixing phase m ; x , the resulting intensity presents os- 
cillations with an amplitude equal to 2 |/ h | on a back- 
ground of value ~ 2|A| 2 iV CT ii nc . 

In conclusion, we have proposed two possible ways 
of detecting a long-range pairing order in a degenerate 
Fermi gas by measuring the coherence function of the 
pairs via matter- wave interferometric techniques. This 
has the advantage over other techniques of directly mea- 
suring the order parameter without relying on a micro- 
scopic description of the many-body state, so that it ap- 
plies in an unambiguous way both in the weakly interact- 
ing and the strongly interacting regimes. More generally, 
the proposed scheme is an application of quantum optics 
techniques to Fermi fields, a line of research expected to 
open new possibilities in the experimental manipulation 
and characterization of fermionic systems. 
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